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Abstract 

In the context of the recent (79.5 kTy) Super-Kamiokande atmospheric neu- 
trino data, we concisely review the status of muonic-tauonic flavor oscillations 
and of the subdominant electron or sterile neutrino mixing, in schemes with 
three or four families and one dominant mass scale. In the three-family case, 
where we include the full CHOOZ spectral data, we also show, through a 
specific example, that "maximal" violations of the one-dominant mass scale 
approximation are not ruled out yet. 

PACS: 14.60.Pq, 13.15.-Fg, 95.85.Ry 
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I. INTRODUCTION 



The Super-Kamiokande (SK) Collaboration has recently presented an updated set of 
atmospheric neutrino data for a detector exposure of 79.5 kTy [|l|. The corresponding 
statistics is about twice as large as compared to the one considered in our earlier published 
analysis of three-flavor mixing p[ and of two-flavor mixing with nonstandard dynamics 
0, and is about 12% larger than in our previous analysis of four-family schemes in 0]. In 
addition, the CHOOZ collaboration presented in final spectral data, which we now include 
in 3z/ analyses 0, improving the accuracy of our previous results p| based on the CHOOZ 
total rate 

Therefore, we think it useful to present a concise update of our 2z/, 3z/, and 4z/ oscillation 
studies, performed under the hypothesis of one mass-scale dominance, so as to elucidate the 
current status and implications of subdominant (electron or sterile) neutrino mixing in SK. 
Finally, we also discuss a specific 3i> example which maximally violates the assumption of 
one-dominant mass scale. 

II. TWO NEUTRINOS 

The standard 2z/ case of z/^j ^ z/^ oscillations can be parametrized by one squared mass 
difference between two states (z^i, z/2), 

= \ml — mil , (1) 
and by one mixing angle ip, describing the flavor content of z/2, 

\{l^2Wp,)\ = , (2) 

|(z/2|z/^)| = , (3) 

where s = sin and c = cos. 

Figure 1 shows our 2i/ best fit to the latest SK data, reached at (m^, sin^ 2ip) = 
(3 X 10^^ eV^,0.97), and corresponding to xLin — 38.5 for 55 — 2 degrees of freedom — a 
very good fit.Q The SK collaboration finds the minimum at a slightly different point, 
(m^, sin^ 2'?/')sK = (2.5 x 10~^ eV^, 1.0) ||l|]. However, the difference is not statistically sig- 
nificant, since the function turns out to be rather flat around the minimum, and the 
"distance" between our best-fit point and the SK one is only about one unit in Ax^. Con- 
cerning the 2z/ bounds on {im?,il)) from parameter estimation, they will be discussed later 
as limits of 3z/ and cases. 

The striking evidence in favor of standard Vr oscillations (Fig. 1) strongly con- 

strains nonstandard explanations. By using the approach described in 0], we parametrize a 



^The overall reduction of Xmin '^ith respect to the 21/ (subcase) analysis presented in Q is mainly 
due to a better agreement of the latest SK electron distributions with their no-oscillation expecta- 
tions (especially in the multi-GeV sample). 
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wide class of scenarios involving nonstandard dynamics through three (free) parameters: an 
oscillation amplitude a, an overall phase factor j3, and an energy exponent n, the standard 
mass- mixing dynamics being recovered for n = —1. 

Figure 2 shows the results of such a three-parameter fit in terms of the projection x^i''^)- 
The corresponding bounds on n give n = — 1.03±0.31 at 90% C.L. (A^^ = 6.25 for three free 
parameters, A^df = 3), in perfect agreement with the standard case. Such results strengthen 
our previous bounds obtained with smaller (45 kTy) SK statistics {n = 0.9 ± 0.4 and 
definitely exclude nonstandard dynamics with integer n 7^ — 1 in the z/^ z/,- channel. 

III. THREE NEUTRINOS 

Oscillations of three neutrinos (z/i, z/2, z/3), under the hypothesis of one mass scale dom- 
inance for atmospheric i/'s (equivalent to set ml ~ m^), are characterized by one squared 
mass difference, 

rri^ = ml- m\^ , (4) 

(the cases > and m? < being physically different [Q), and by two mixing angles 
{ip = 6*23 £ [0, vr/2] and = 6^13 G [0, 7r/2]), describing the flavor content of the state 1/3, 

|(z/3|z/e)| = s<^ , (5) 
\W3W,,)\ = c<ps^ , (6) 

\{i^3Wt) \ = c^c^ , (7) 

the pure z/^ ^ Ur case being recovered for = (see and references therein). 

In order to show the bounds in the mixing parameter space, we used in [0] a triangular 
representation |^,rD|, which is basically a linear mapping in the variables (sin^ sin^ 0). 



Such a representation has the advantage of embedding unitarity by construction, but has 
the disadvantage of not showing in detail the phenomenologically interesting case of small 
(J). We use here the alternative representation in terms of (tan^ ip, tan^ 0) in logarithmic 
scale (also introduced in [^|lOl), which expands the small region while preserving octant 
symmetry (when applicable). 

Using SK data only, and assuming > 0, we find the best fit (Xmin — 38.1) at 
(m^jtan^^/'jtan^^) = (3 x 10"^ eV^, 0.9, 0.01). The slight deviation of the best fit mix- 
ing from the pure i/^ ^ Ur maximal mixing [(tan^ ^Z', tan^ 0) = (1,0)], although intriguing^, 
is — unfortunately — not statistically significant (Ax^ ^ 1). This also implies that there is no 
significant indication for possible matter effects related to z/g mixing (tan^ > 0) in the SK 
data. 

Figure 3 shows the 3z/ volume allowed at 90% and 99% C.L. (Ax^ = 6.25 and 11.36 for 
Ndp = 3, respectively) in the (m^, tan^ ■0) tan^ 0) parameter space, through its projections 
onto the coordinate planes. The upper limit on tan^ (< 0.35 at 90% C.L.) improves the 



^There is no reason to have exactly (tan^ ^,taii^ 0) = (1, 0). 
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one found in (tan^ < 1 at 90% C.L. from 33 kTy SK data), showing the steady progress 
of SK in confirming dominant Uj. mixing and in constraining additional mixing. The 

90% C.L. range for is (1.6-7.2) x 10~^ eV^. The bounds on tan^ ip in Fig. 3 are octant- 
symmetric only in the 2z/ limit tan^ (as they should), and show a slight preference 
for ip in the second octant when tan^ > 0. Correspondingly, slightly higher values of 
are preferred. The (weak) positive correlation between tan^ and tan^ or m^, however, is 
largely suppressed by the inclusion of CHOOZ data, as we now discuss. 

As described in [0, we can now include the CHOOZ reactor spectral data (14 bins minus 
one adjustable normalization parameter) through a statistics reproducing the bounds of 
the so-called "CHOOZ analysis A" |^. This improvement provides more accurate bounds 
in the region of interest for atmospheric neutrinos. Our best fit to SK+CHOOZ data 
(Xmin = 45.7) is reached at (m^, tan^ ^, tan^ 0) = (3 x 10"^ eV^, 0.75, 0.003). Once again, 
the small deviation of the best fit mixing from (tan^ ■0, tan^ 0) = (1,0) is not statistically 
significant (Ax^ ^ 1). 

Figure 4 shows the projections of the (m^, tan^ -0, tan^ 0) volume allowed by 
SK+CHOOZ. By comparing Fig. 3 with Fig. 4, the tremendous impact of CHOOZ on 
Uf. mixing bounds becomes evident (one order of magnitude difference in the upper bound 
on tan^ 0). As expected, at the small values of tan^ allowed by the fit in Fig. 4, both the 
octant-asymmetry in ip and the upper limit on are reduced, and the 90% C.L. range for 

becomes (1.6-5.3) x 10~^ eV^. 

We have also repeated the fit for the case < (not shown), corresponding to a state 
1/3 lighter than i>i^2- For negative m^, we get somewhat weaker bounds on tan^ (^ 0.5 at 
90% C.L.) for the fit to SK data only, while the fit to SK+CHOOZ data gives results almost 
identical to those in Fig. 4. This fact shows that, unfortunately, current atmospheric+reactor 
data are basically unable to discriminate the sign of in 3i> scenarios, as it was the case 
for pre-SK and pre-CHOOZ data |[. 

Finally, Fig. 5 shows the SK zenith distributions computed for three representative 
cases at tan^ = 0.025 (allowed at 90% C.L. by SK+CHOOZ) and for ip both maximal 
[tan'^ijj = 1) and nonmaximal (tan^'0 = 1/2 and 2). Within statistical errors, the Su zenith 
distributions in Fig. 5 are hardly distinguishable from the 2z/ one in Fig. 1 (even more so for 

< 0, not shown), the differences being at most ~ 1.5a in a few bins. Therefore, there is 
little hope to unambiguously discover 7^ (i.e., Ue mixing) from SK atmospheric data in 
the near future. 

IV. FOUR NEUTRINOS 

We consider the Au (3 active + 1 sterile) scenario described in [Q, characterized by a 
2+2 mass spectrum with well-separated atmospheric and solar doublets. We also make the 
simplifying assumption ^ that the atmospheric doublet (1^3,2/4) is almost decoupled from 
Uf,, and that the solar doublet {ui, z/2) is almost decoupled from u^. The dominant mass scale 
for atmospheric neutrinos is then 

= ml — ml , (8) 
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and two mixing angles {ip = 623 G [0,7r/2] and ^ = ^13 G [0,7r/2]) are sufficient to describe 
the flavor contents of the state u^, 

|(l/4|z/e)| ^0 , (9) 

\{j^4Wt^)\ = , (10) 

\{l^4\l^r)\ = C^C^ , (11) 

|(z/4|z/,)| = C^Sg , (12) 

the pure z/^ Uj. and z/^ ^ z/^ cases being recovered for = and = 1, respectively. In 
such scenario, the cases > and < are not physically different (being equivalent 
under octant inversion, ip —* 7t/2 — ip 0]), and we take > 0. 

Figure 6 shows the bounds obtained by a fit to SK data at 90% and 99% C.L. (A^df = 3) 
in the parameter space (m^, tan^ ■?/', tan^ ^). The best fit point (Xmm = 38.1) is reached at 
(m^, tan^ '?/', tan^ (^) = (3 x 10~^ eV^, 0.76, 0.1) but, once again, the preferred mixing differs 
from (tan^ tan^ 0) = (1,0) by less than one unit in Ax^. Notice that in the (m^, tan^ ■?/') 
plane, the projected bounds of Fig. 6 are very similar to the 3z/ ones in Fig. 4, implying 
that the current limits on the "2z/" subset of parameters (m^, tan^ ■?/') are rather stable even 
by making allowance for additional Ug or i/^ mixing. 

Concerning the (tan^ ^z^, tan^ ^) plane in Fig. 6, the upper bounds on tan^^ indicate 
that pure ^ Vg oscillations (tan^ 00) are disfavored as compared with pure z/^ ^ z/^ 
oscillations (tan^ ^ — 0), in agreement with [^,0, although large Vs mixing is not excluded 



yet. In particular, the current upper limit from Fig. 6 (tan^ ^ ^ 4 at 90% C.L.) is even slightly 
weaker than the one we found with smaller statistics in |^ (tan^^ ^ 2). The reason can 
be traced to a peculiar feature of the latest UP/i data, namely, the flatness of the muon 
suppression pattern in the four UP/x bins at cos^^ G [—0.7, —0.4], as described in Fig. 7. 

Figure 7 shows the SK zenith distributions for three representative 4z/ cases with sizable 
Vs mixing (tan^,^ = !)• As it is well known (and evident from a comparison of Fig. 7 with 
Fig. 1), additional Vg mixing for atmospheric neutrinos tends to reduce the muon suppression 
and, in particular, tends to flatten the normalized UP/i distribution. Although the SK UP/x 
data do prefer a mean positive slope rather than a flat suppression, the four bins in the 
zenith range cos9 G [-0.7,-0.4] happen to favor a locally flat distribution. This current 
feature might be just a statistical fluctuation but, at present, it plays some role in global 
fits, where it tends to weaken the rejection of "flat" distributions (i.e., of sizable Vg mixing), 
as compared with previous UP/z data ||^. 

The SK Collaboration has also presented additional (preliminary) indication in favor of 
z/^ ^ mixing coming from statistical z/^ appearance in selected event samples [|T|. It 
seems possible to isolate an excess of about 100 ± 50 r-like events, to be compared with 
standard z/^ ^ Vr expectations of ~ 100 |]l|. Taken at face value, such numbers imply an 
additional ~ 2a evidence (Ax^ ~ 4) in favor of pure z/^ (s| = 0) as compared with 

pure z/^ ^ Vg (^l = 1). We have then roughly parametrized the SK "tau appearance" signal 
by adding a penalty function Ax^ = 4s| in the 4z/ fit. We get an "improved" upper bound 
tan^^ ^ 1.5 at 90% C.L. (A^df = 3), to be compared with tan^^ ^ 4 in Fig. 6 (without 
penalty function). On the one hand, this seems to indicate that there is certainly room 
to refine current bounds on additional Vg mixing in SK; on the other hand, our analysis 
shows that large Vg mixing (e.g., a fifty-fifty admixture of z/^ and Vg at tan^^ ~ 1) is not yet 
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excluded at present. Therefore, compatibility with complementary solar neutrino bounds 
on tan^^ is still possible, as it was the case for previous SK data [Q. 



V. TWO, THREE, AND FOUR NEUTRINO SUMMARY 

The bounds on dominant and subdominant mixing found in the previous sections can 
be conveniently summarized in one single plot, as shown in Fig. 8. The left panel shows 
the bounds on the (m^, tan^ ■?/') parameters for pure z/^ Vj. mixing (equivalent to our 
scheme for tan^ = or to our 4z/ scheme for tan^ ^ = 0). As discussed before, such bounds 
are not significantly altered by additional z/g mixing (tan^ > 0) or by additional Vg mixing 
(tan^^ > 0), and thus they hold also in the global 3i/ and 4i/ fits with good accuracy. For 
such reason, the 2v bounds in the left panel are formally obtained for iVop = 3, so as to 
match those in the middle and right panels. The middle panel shows the 3z^ bounds on 
additional v^. mixing (tan^ > 0), with and without CHOOZ. Finally, the left panel shows 
the bounds on additional Vs mixing (tan^^ > 0). Such synthetic figure represents the 
main result of our analysis. 

VI. TWO MASS SCALES 

The analyses in the previous sections are based on the assumption that atmospheric 
neutrino oscillations are driven by only one mass scale {m?). This is not necessarily the 
case, especially if one takes the "solar" squared mass difference in the upper range allowed 
by the data (see, e.g., [Q), provided that one accepts an averaged or quasiaveraged solar 
neutrino survival probability 0. Concerning atmospheric neutrinos, two mass scales in the 
range ~ 10~^ eV^ were shown to provide acceptable fits to previous SK+CHOOZ data | |15[| . 
Here we show, through a specific example, that such possibility is not yet excluded by the 
latest data. 

Let us consider the specific 3z/ case shown in Fig. 9, characterized by two equal 
squared mass differences {Am\2 = ^fn^i = 0.7 x 10~^) eV^, and by mixing angles^ 
(tan^ 0, tan^ '0, tan^ = (0, 1,2), giving the flavor composition of mass eigenstates shown 
in the same figure. The spectrum in Fig. 9 might be called "democratic" , since it maximally 
violates the usual "hierarchical" approximation {Am2i <^ Am^g). Notice that z/g oscillations 
are driven only by Am^i, which is purposely chosen just below the current CHOOZ bounds 
H. The solar neutrino survival probability is then P(z/e ^ z/e) ^ 1- |sin2 2a; ~ 5/9 ~ 1/2, 
up to small (quasiaveraged 0) corrections. 

Figure 10 shows the SK zenith distributions computed for the democratic scenario in 
Fig. 9, and corresponding to = 50.9 for the fit to SK data only [x^ = 61.2 if CHOOZ 
data are also included). Although such value is significantly higher than in the best-fit 2u 
case of Fig. 1, it is still acceptable from the point of view of goodness of fit. All in all, 
the curves in Fig. 10 provide a globally acceptable "fit-by-eye", with moderate departures 



^In standard notation, = ^13, ip = 623, and to = 612- 
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from the data along the horizontal direction (cos 9^0) for the MG/x and US/i samples. 
Therefore, if one accepts an almost constant (~ 1/2) suppression as explanation of the solar 
neutrino deficit, scenarios with two comparable mass scales appear to represent a viable 
possibihty in the current atmospheric neutrino phenomenology. Needless to say, the most 
general analysis of such cases (depending on all the Su mass-mixing parameters) would be 
rather intricated, and is postponed to a future work. 

VII. CONCLUSIONS 

In the context of the latest (79.5 kTy) SK atmospheric u data, we have concisely reviewed 
the status of dominant ^ oscillations (including the case of nonstandard dynamics. 
Fig. 2) and of subdominant and mixing (in and 4i/ schemes, respectively). In 
the ?>v case we have applied an improved CHOOZ analysis. The main 2z/, 3i/, and Ai/ 
results are discussed separately and then summarized in Fig. 8. Finally, we have shown 
(through a specific example) that current atmospheric neutrino data are also compatible 
with oscillations driven by two comparable mass scales. 
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Fig. 1. Super-Kamiokande zenith distributions (79.5 kTy Q) used in the analysis, normahzed 
to no-oscihation expectations. The data set includes sub-GeV electrons (SGe, 10 bins), multi- 
GeV electrons (MGe, 10 bins), sub-GeV muons (SG/i, 10 bins), multi-GeV muons (MG^, 10 bins), 
upward stopping muons (US^, 5 bins), and upward through-going muons (UT/i, 10 bins), for a total 
of 55 data points. The error bars are statistical only (itlcr); systematic (correlated) uncertainties 
are treated as in p. The solid line is our best fit for 2v oscillations {x^ = 38.5). 
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Bounds on nonstondord dynamics 
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Fig. 2. Dependence of the function on the neutrino energy exponent n, assuming an oscil- 
lation phase proportional to E^, with unconstrained factors for the overall phase and amplitude. 
The only integer n compatible with the SK data is n = — 1, corresponding to standard mass-mixing 
dynamics. 
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Fig. 3. Projections of the regions allowed in the parameter space {m , tan^ V, tan^ (p) at 90% 
and 99% C.L. (Ax^ = 6.25 and 11.36 for A^df = 3) onto the coordinate planes. The fit includes 
SK data only (79.5 kTy). The pure 2u case of <->• 1/^ oscillations is recovered for tan^ — >■ 0. 
Nonzero values of (f) parametrize Ve mixing. 
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Fig. 4. As in Fig. 3, but including final CHOOZ positron spectra Q (14 data points minus one 
adjustable normalization factor). 
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Fig. 5. Zenith distributions for three representative 3z/ cases with tan^ cj) = 2.5 x 10 ^, ahowed 
at 90% C.L. by SK+CHOOZ. Notice the distortion of the MGe distribution. Such distortion would 
be somewhat smaller for negative = — 3 x 10~^ eV^ (not shown). 
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Fig. 6. Projections of the regions allowed in the 4i/ parameter space (m^, tan^ i/;, tan^ ^) at 90% 
and 99% C.L. (Ax^ = 6.25 and 11.36 for A^df = 3) onto the coordinate planes. The fit includes 
SK data only (79.5 kTy). The pure case (tan^^ 0) is clearly preferred over the pure 

case (tan^^ — > co). However, sizable Vs mixing [tan^ ~ 0(1)] is allowed in addition to 
t'it ^ oscillations. 
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Fig. 7. Zenith distributions for three representative Au cases with sizable Ug mixing (tan^^ = 1), 
ahowed at 90% C.L. by SK data. Notice the reduced suppression in the muon samples, as compared 
with pure v-j- mixing in Fig. 1. 
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Fig. 8. Summary of Sf, and 4z/ bounds at 90% and 99% C.L. on the mass-mixing parameters 
from SK data. Left panel: Bounds on (m^,tan^V) for pure Uj- mixing (i.e., for tan^ ^ = 

= tan^^). Middle panel: Bounds on additional mixing (parametrized by tan^ (6 > 0) in 
scenarios, both without CHOOZ (sec also Fig. 3) and with CHOOZ (see also Fig. 4). Right panel: 
Bounds on additional Vg mixing (parametrized by tan^ ^ > 0) in 4z/ scenarios (see also Fig. 6) . All 
the bounds are derived for Ndf = 3, including those in the left panel. 
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'Democratic" scenario 
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Fig. 9. A "democratic" (i.e., nonhierarchical) scenario with two equal squared mass 
differences — ^"^-ii — ^•'^ ^ 10~^ eV^) and with {i'(.,v^,Vt) flavor content as follows: 

(1/3,1/3,1/3) for vi, (1/6,1/6,1/6) for 1^2, and (0,1/2,1/2) for 1/3. Notice that i^e disappearance 
is driven only by Am2i (just below the CHOOZ sensitivity [^), and gives P(z^e — ^ i^e) — 5/9 for 
solar neutrinos (up to small quasiaveraged oscillation corrections). 
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Fig. 10. Zenith distributions for the democratic scenario in Fig. 9, giving 
to SK data (x^ = 61.2 for the fit to SK+CHOOZ). 
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